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We study quantum transport in Dirac materials with a single fermionic Dirac cone (strong topolog- 
ical insulators and graphene in the absence of intervalley coupling) in the presence of non-Gaussian 
long-range disorder. We show, by directly calculating numerically the conductance fluctuations, that 
in the limit of very large system size and disorder strength, quantum transport becomes universal. 
However, a systematic deviation away from universality is obtained for realistic system parame- 
ters. By comparing our results to existing experimental data on 1// noise, we suggest that many 
of the graphene samples studied to date are in a non-universal crossover regime of conductance 
fluctuations. 



With the successful isolation of graphene pQ, and the 
discovery of topological insulators (TI), the study of 
Dirac materials has come to the fore. In a Dirac material, 
the conduction and valence bands touch at an isolated 
set of points, called Dirac points (DP) (or valleys). The 
energy spectrum close to the DPs is linear and the low 
energy properties of these materials are well described 
by a set of massless 2D Dirac fermions. Graphene is in- 
trinsically a 2D material [2], while in the case of a 3D 
topological insulator the Dirac fermions are realized as 
surface states whose presence is demanded by the non 
trivial topology of the bulk energy bands [3]. In this 
Letter, we report an exact numerical theoretical study 
of quantum conductance fluctuations (CF) in the pres- 
ence of non-Gaussian long-range disorder in 2D Dirac 
materials, which should apply directly to both graphene 
(in the absence of intervalley scattering) and TI surface 
states. Specifically by (non-Gaussian) long-range disor- 
der we refer to disorder whose spatial correlation function 
decays algebraically. In conventional parabolic band 2D 
electron liquids, taking into account screening, the dis- 
order potential is short-range even when due to charge 
impurities. In Dirac materials, in contrast, due to the 
massless nature of the fermionic excitations, the disorder 
potential due to charge impurities retains its long-range 
nature even when screening is taken into account; in par- 
ticular the spatial correlation function decays as 1/r 3 at 
large distances. As a consequence our work is the first to 
address the very fundamental and broad question of how 
the long-range nature of the disorder affects the quantum 
transport properties of 2D electron systems. 

In graphene the dominant source of disorder is believed 
to be from remote charge impurities [2] (though other 
scattering mechanisms might also be of relevance in cer- 
tain cases [4]). This type of disorder does not couple 
the two valleys of graphene and to a good approximation 
graphene can be described by a single Dirac fermion. At 



the same time, a strong TI is characterized by an odd 
number of Dirac cones in the surface spectrum (as op- 
posed to a weak TI which has an even number). The 
dominant physics is captured by the case of a single Dirac 
fermion which is experimentally realized in e. g. Bi2Sc3. 
In both cases the Hamiltonian is invariant under a time 
reversal symmetry T with T 2 = —1, placing these sys- 
tems in the symplectic (All) symmetry class [5]. Below 
we focus on the single Dirac fermion case. 

In the clean case, the density of states goes to zero and 
transport is by quantum tunneling of evanescent modes 
rather than diffusion of propagating modes. This regime 
is referred to as pseudo- diffusive [6j. Disorder drives the 
system away from the pseudo-diffusive regime and into 
a symplectic metal phase characterized by weak anti- 
localization [7 (enhanced conductivity due to destructive 
interference of time reversal symmetric paths [8 ). 

These conclusions were reached by studying the effects 
of short range disorder. However in a Dirac material the 
disorder potential created by charge impurities retains its 
long range character even after screening has been taken 
into account [5]. This leads to strong carrier density in- 
homogeneities and, close to the DP it induces puddles of 
electrons and holes. In this regime, conventional analyti- 
cal methods fail and one needs to rely on numerical sim- 
ulations. This puddle formation due to random charged 
impurities in the environment has been experimentally 
observed both in graphene [TUj and TIs [IT] . 

In this paper we study the problem of quantum trans- 
port in the presence of long range disorder both in the 
puddle regime and away from the DP at large carrier den- 
sity, n. We study both quantum corrections to the aver- 
age conductivity (a) [weak anti-localization (WAL)] and 
the CF arising from changes in interference patterns as 
external parameters are varied. In conventional semicon- 
ductors and metals the amplitude of these fluctuations 
is a universal number giving rise to universal conduc- 
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tance fluctuations (UCF) . In terms of the scaling function 
j3(a) = dlna/dlnL, where L is the sample length, the 
WAL is the first order term in the expansion in \ ja. It 
happens that the second and third order terms are iden- 
tically zero in the symplectic class [T2] and therfore the 
WAL can be obtained already for rather small values of 
cr. In contrast, higher order terms in the UCF are ex- 
pected to be nonzero making the numerical calculation 
of UCF more sensitive to finite size effects. Moreover, 
in general, the calculation of CF requires a much larger 
number of disorder realizations than the calculation of 
ex. These facts make the numerical calculation of CF 
in the symplectic class much more challenging than the 
calculation of a and WAL corrections. Earlier numeri- 
cal studies [T3] of UCF in graphene using tight binding 
models with short range disorder found results consistent 
with the UCF theory [T3], however, even for short-range 
disorder, no systematic study of density dependence or 
deviations from universality was attempted. 

We demonstrate that just as in the case of short range 
Gaussian disorder, long range disorder drives the system 
into the symplectic metal phase. However, the crossover 
regime before the universal quantum transport sets in is 
considerably larger than in the case of short range dis- 
order. This gives rise to CF with an amplitude which is 
smaller than the UCF close to the DP but considerably 
larger at high n which, as we argue, can lead to an in- 
triguing non-monotonic behavior of the 1 // noise, which 
may have already been observed in graphene. We em- 
phasize that our study is performed at zero temperature, 
in the absence of inelastic scattering, conditions in which 
the phase coherence length is infinite. Our analysis 
is valid as long as Z p h > L. The main goal of our work 
is to study the effect of the long-range nature of the dis- 
order on the universal properties of quantum transport, 
for this reason we must have l p h > L otherwise sample- 
dependent effects due to finite i p h would mask the sought 
after effects. Moreover, any quantum transport measure- 
ment strives to be in this regime that appears also to be 
the relevant one for current experiments on mesoscopic 
samples of graphene. 

To obtain the conductance G of a single Dirac fermion 
we solve, using the transfer matrix method described in 
Ref. [7J, the scattering problem defined by the Hamil- 
tonian: H = vpp ■ cr + Vo(r) — /i. Here vf is the 
Fermi velocity, cr the vector formed by the Pauli ma- 
trices (<T x ,a y ) (in spin space for TIs and in sublattice 
space for graphene) , [i the chemical potential which con- 
trols n, and Vp is the disorder potential. We take a 
sample of length L and width W with periodic bound- 
ary conditions in the transverse direction. The solution 
of the scattering problem returns the transmission am- 
plitudes t that in turn give us the two-terminal conduc- 
tance G = R^ 1 = (g s g v e 2 /h)Tr(tH), R being the re- 
sistance and g s and g v the spin and valley degeneracies 
respectively (g s g v = 4 for graphene and g s g v = 1 for the 



surface of a TI). In the remainder for concreteness we 
set g s g v = 4. The conductivity is obtained via the rela- 
tion a = {WdR/dL)- 1 which minimizes the contribution 
of contact resistance [T3]. The CF are in turn obtained 
as the variance ((SG) 2 ) = (G 2 ) - (G) 2 , where the angle 
brackets denote the average over disorder realizations. 

The scattering potential Vd is induced by the remote 
impurity charges. In the absence of screening, in momen- 
tum space, Vo(q) = (e 2 /e n )A(q)e~ qd /q, where eo is the 
background dielectric constant, d the average distance 
of the impurities from the graphene layer, and A(q) are 
random numbers with Gaussian distribution such that 
(A) = and (A 2 ) = n imp , with ni mp the impurity density. 
The bare disorder is renormalized by the interactions re- 
placing the dielectric constant eo by a dielectric function 
e —> e(q). To calculate e(q) in the presence of charge im- 
purities we use the Thomas Fermi Dirac Theory (TFD) in 
which both Hartree and exchange-correlation terms are 
taken into account [THJ [T7J . The TFD approach is valid 
when the inequality |Vn|/|n| < Uf is satisfied. Close to 
the DP, due to the disorder, the root mean square of n 
(n( rms )) is much larger than its average. In this situa- 
tion n' rms ) should be used instead of n in the inequality 
above to verify the validity of the TFD theory [THl HH] • 
Because n( rms > « n imp [H E5] the TFD theory is valid 
as long as n; mp is not too small (> 10 n cm~ 2 ) |2"U] . 
The TFD gives results for the carrier density profile n(r) 
and e(q) that compare well with both density functional 
theory (DFT) results [ST] and experiments [TU] and is the 
only method available, for large samples and ensembles, 
that takes into account nonlinear screening effects that 
dominate close to the DP. 

For purposes of comparison it is useful to consider dif- 
ferent models of screening. In particular one can turn off 
exchange-correlation terms in the TFD. We denote the 
corresponding dielectric function by e nxc , the full TFD di- 
electric function being ctfd ■ The cases of doping depen- 
dent but uniform screening e _1 — > e7 d = q/(q + 4r s kp), 
where r s — e 2 /(cqJivf), and constant screening e _1 — > 
e cs = <?/(<? + Is) w iU also be considered. Finally, we also 
compare with the case of Gaussian correlated disorder de- 
fined by (V D {v)V D {v')) = /f /(2< 2 ) exp(-|r- rf /2£ 2 ). 
K Q is a dimensionless measure of the disorder strength 
and £ the correlation length. This is the type of disorder 
correlation considered in most prior numerics, and in an- 
alytical considerations (where the limit £ — > is routinely 
taken) . 

The analytical theory of UCF requires diffusion, i.e. 
L ^> £, where £ is the mean free path. In this limit 
one can perform a controlled diagrammatic calculation of 
transport properties. Adapting the results of the short- 
range calculation in [14] to our geometry, this approach 
predicts a universal (independent of microscopic param- 



3 




n(10 12 cm" 2 ) n. (10 12 cm 2 ) 

* ' imp * ' 

FIG. 1. (Color online), (a) (a) as a function of n for ni mp = 
(3,5,9) x 10 12 cm" 2 from top to bottom, (b) (a) at the DP 
as a function of ni mp for the case of Coulomb disorder. The 
inset shows (a) as a function of sample length L = x for W = 
160 nm, and n = (3,2,1) x 10 12 cm -2 from top to bottom. 
The dashed lines show da/dlnL — 4e 2 /Trh. 




FIG. 2. (Color online), (a) ((SG) 2 } as function of n for the 
case of long-range disorder (L = 160 nm) and different values 
of ni mp ; from top to bottom (for n = 2 x 10 12 cm -2 ): rii mp = 
(0.682, 1, 3, 5, 9, 16) x 10 12 cm~ 2 . (b) ((SG) 2 ) as a function of n 
for different types of disorder potential (see text), (c) {(SG) 2 } 
as a function of kp^ for Gaussian disorder, (d) ((SG) 2 ) at the 
DP as a function of ni mp for the case of long-range disorder 
(L = 160 nm). In the inset {(SG) 2 } as a function L/W; the 
top line shows the results for ni mp = 5 x 10 12 cm~ 2 and 
n = 3x 10 12 cm -2 ; the solid line below, and the bottom one, 
show the results for the same n lulp and n — 10 12 cm -2 and 
n = respectively. Between these two lines there is the line 
showing the results for the case of Gaussian disorder with 
A"o = 1 and = 0.8. The dashed lines in all the panels 
show the universal value of Eq. (TTj> 



eters, such as n- lmp , Kq and £) value of the CF given by 
Fig. [lja) shows the dependence of (a) on the doping n 



for three different values of n; mp . (a) approaches linear 
dependence on n at large densities in qualitative agree- 
ment with previous theoretical results [18j [22] and ex- 
periments [53]. At the DP, the conductivity (cr m j n ) is 
enhanced by disorder, see Fig. [TJb), in contrast to what 
is obtained in Boltzmann type of theories extended to the 
DP, but in agreement with simulations with short range 
disorder. For large system sizes and disorder strengths 
a increases logarithmically with system size, as shown 
in the inset of Fig. [ijb), consistent with weak anti- 
localization which predicts a ~ (g s g v /Tt)bxL |24j . 

Fig. ga) shows ((SG) 2 ) as a function of n for several 
values of rii mp . The dashed lines in Fig. [2] show the uni- 
versal value of Eq. . We observe that close to the DP 
the CF tend to be smaller than the UCF value, while at 
large n they are larger. As n- lmp is increased, the CF ap- 
proach the universal value with the set of curves crossing 
roughly at the UCF value. Similar trends are observed as 
a function of system size. This suggests that in the ther- 
modynamic limit the CF in the presence of long range 
disorder, approach the universal value independent of n. 
In this limit the diffusive symplectic metal is an accu- 
rate description of the system. In contrast, percolation 
predicts CF larger than UCF at the DP [25] . 

To understand better the deviation from universal- 
ity and to quantify the behavior observed in Fig. [2] 
we estimate £ using the Boltzmann theory, which gives 
£ = f{r s ,kwd)^/n/ni m p> 0, where / is a function that 
for d < 1 nm depends very weakly on kpd [22 [2S]. For 
d = 1 nm and r s — 0.8, / ~ 5.66 almost independent 
of hp- In Fig. [2] we represent with open symbols those 
data points for which £ > 0.25L, revealing that the ob- 
served oscillations at large n and small n; mp are ballistic 
effects. It also suggests that the departure of the CF 
from the UCF value at high n is due to the fact that we 
are in the ballistic to diffusive crossover regime. Simi- 
larly, close to the DP, the deviations are attributed to a 
pseudo-diffusive to diffusive crossover (the estimates of £ 
are not as reliable in this regime). 

To better understand the role played by the disorder- 
induced density inhomogeneities we plot in Fig. [2jb) the 
CF for the different variants of screening discussed above. 
We observe that neglecting exchange-correlation terms 
in the TFD increases the deviation of the CF from the 
universal value close to the DP. This is due to the fact 
that in the absence of exchange-correlation the amplitude 
of the density fluctuations are increased [16]. In con- 
trast, assuming uniform screening, either doping depen- 
dent or constant, gives results that are closer to the UCF 
value. The results of Fig. |2jb) clearly show that the pres- 
ence of strong carrier density inhomogeneities increases 
the range of n and n; m p for which the transport is in 
the crossover regime for which the CF differ from the 
universal value. 

In Fig. [2jc) we show the CF for the gaussian corre- 
lated potential as a function of kp£, — (^nn/gsgv) 1 ^ 2 . 
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In this case from the Boltzmann theory we have I = 
2y/2ir{; 3 kp/K Q . As in the case of charge impurities we 
find that in the pseudo-diffusive regime ((5G) 2 ) is smaller 
than the UCF value, and larger in the ballistic high den- 
sity regime, but approaches the UCF value as the im- 
purity strength Kq or system size L are increased. In 
contrast to the long range case, there is a large range of 
intermediate values of n where the CF quickly goes to 
the UCF value. Thus, the crossover regime is strongly 
suppressed for short-range disorder compared with long- 
range disorder. In the presence of spatial correlations 
among charge impurities the long-range nature of the 
disorder is suppressed compared to the case of uncorre- 
cted impurities [27] and so we predict that the crossover 
regime will also be reduced. 

The trend that emerges from Fig. [2^ a) is that as 7ii mp 
increases the CF approach the universal value. This is 
further demonstrated in Fig. [27d) which shows the value 



(a)' 



of the CF at the DP as a function of ni mp for the case 
of charge impurities. We see that for very large nj mp , 
for which t/L « 1, the CF saturate to the universal 
value (dashed line). In the inset to Fig. [5Jd) we compare 
the aspect ratio dependence of the CF to the analytical 
expression ([!]). Both in the case of long range and Gaus- 
sian disorder the curves agree well within the universal 
regime. In the crossover regime the CF follow the same 
trend but with an amplitude that differs from the UCF 
prediction. 

Taken together these results suggest that even in the 
presence of long-range disorder, the system is eventually 
(large enough L and/or disorder strength) driven into the 
universal symplectic metal phase and that, however, the 
crossover regime in which neither ballistic or universal 
diffusive physics is applicable is very large in the presence 
of long-range Coulomb disorder and puddles. 

To connect these results with experiments we consider 
the relevance of the CF results to the 1// noise. In 2D 
(and ID) the displacement of a single defect can cause a 
change in G of the order oie 2 /h [25] . One of the probable 
sources of resistance fluctuations is thermally activated 
motion of defects. Assuming that the defects move on a 
time scale t S> Ti no i astic and that the hopping distance is 
uncorrelated with r, for the spectrum of the resistance 
fluctuations Sr(w) = J dt(5R(t)SR(0))e iut we have [28] 



Sr{w) = R 2 



((SG) 3 
(G) 2 



2r 



+ 1 



P(r)dr (2) 



where P(t) is the probability distribution of the time 
that it takes a defect to move. If the motion of defects is 
thermally activated then P(t) is quite broad and Sr(u>) 
will be approximately l/f. We can then calculate the 
dependence with respect to the doping of the strength of 
the l/f noise by calculating the coefficient ((5G) 2 ) / (G) 2 , 
or equivalently ((SG) 2 ) / (a) 2 . 

Fig. Hj shows ((8G) 2 )/(a) 2 as a function of n for both 
charge disorder and Gaussian disorder. In both cases 




FIG. 3. (Color online), (a) {{SG) 2 )/G 2 as function of n for 
the case of long-range disorder (L = W = 160 nm) and differ- 
ent values of ni mp ; from top to bottom (for n = 2xl0 12 cm -2 ): 
n imp = (5,3, 1,0.682) x 10 12 cm~ 2 . (b) {(5G) 2 )/G 2 as a func- 
tion of &;f£ for the case of Gaussian disorder for Kq = 4, top, 
and K = 1, bottom. L = W = 75£. 



we see that at low disorder strengths ((5G) 2 ) / (a) 2 de- 
pends non-monotonically on n. The reason is that in 
the pseudo-diffusive regime ((SG) 2 } grows rapidly with 
n while G grows very slowly, whereas at larger values of 
n the transport is diffusive and ((SG) 2 ) is almost con- 
stant while G grows rapidly with n. In contrast, in 
the symplectic metal regime the ratio ((SG) 2 ) / '(c) 2 de- 
creases monotonically with n. We conclude that the non- 
monotonic dependence of the l/f noise is a generic prop- 
erty of the crossover regime between pseudo-diffusive and 
diffusive regimes. The l/f noise of graphene has been 
measured experimentally |29j and its dependence on n 
has been shown to qualitatively follow the non-monotonic 
behavior shown in Fig. [3j This suggests that the samples 
in these experiments are not fully in the diffusive regime 
but rather in a non-universal crossover regime. 

In conclusion, we have calculated the amplitude of the 
conductance fluctuations in Dirac materials with a sin- 
gle Dirac cone and in the presence of disorder with var- 
ious range and screening properties. Our work is the 
first to consider the effect of the long-range nature of 
disorder on the quantum transport properties of 2D elec- 
tron systems. In the limit that the mean free path is 
much smaller than the system size, the CF approach the 
universal value predicted by diagrammatic calculations. 
Before reaching the universal value the CF are systemati- 
cally smaller than the universal value close to the DP but 
larger away from it. In particular, for system parameters 
realistic to graphene, the CF seems to deviate from the 
universal value leading to a non-monotonic dependence 
of l/f noise, as recently observed in experiments |29j . 
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